基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

Cutting path-dependent machinability of SiCp/Al composite under multi-step ultra-precision diamond cutting

Particle-tool interactions, which govern the synergetic deformation of SiC particle reinforced Al matrix composites under mechanical machining, strongly depend on the geometry of particle position residing on cutting path. In the present work, we investigate the influence of cutting path on the machinability of a SiCp/Al composite in multi-step ultra-precision diamond cutting by combining finite element simulations with experimental observations and characterization. Be consistent with experimentally characterized microstructures, the simulated SiCp/Al composite is considered to be composed of randomly distributed polygonally-shaped SiC particles with a volume fraction of 25vol%. A multi-step cutting strategy with depths of cut ranging from 2 to 10 μm is adopted to achieve an ultimate depth of cut of 10 μm. Intrinsic material parameters and extrinsic cutting conditions utilized in finite element simulations of SiCp/Al cutting are consistent with those used in corresponding experiments. Simulation results reveal different particle-tool interactions and failure modes of SiC particles, as well as their correlations with machining force evolution, residual stress distribution and machined surface topography. A detailed comparison between numerical simulation results and experimental data of multi-step diamond cutting of SiCp/Al composite reveals a substantial impact of the number of cutting steps on particle-tool interactions and machined surface quality. These findings provide guidelines for achieving high surface finish of SiCp/Al composites by ultra-precision diamond cutting.     

柔性支撑控制力矩陀螺动力学建模与控制

针对柔性支撑引起控制力矩陀螺(CMG)框架控制性能下降的问题，建立了柔性支撑条件下控制力矩陀螺的动力学模型，分析了柔性支撑扰动力矩的耦合机理并给出了解析模型.据此，设计了一种积分滑模控制策略，用于提升对CMG框架动力学特性变化和新增扰动力矩的鲁棒性.仿真结果表明，该控制策略相比传统的PID控制具有更好的动态过程和稳态性能，提升了柔性支撑下控制力矩陀螺框架速度的闭环控制性能.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

一种基于VPL 的元器件快速建模方法

为了更好地落实航空电子模块的可制造性设计审查工作,结合航空电子模块产品特点和技术要求,建设了基于Valor软件的PCB/PCBA设计与检测自动化分析系统.元器件模型库作为该系统的必要组成部分,其覆盖范围和数据精度对审查结果至关重要.本文将从元器件模型库建设必要性、建设过程和建模工艺要求三个方面进行阐述,并提出一种基于Valor VPL的元器件快速建模方法,应用结果表明,该方法能够高效构建满足系统审查精度的元器件模型,使大规模元器件模型库的建设成为可能.     

Development and performance evaluation of tropospheric scintillation model on Ku-band satellite link over Akure,Nigeria

Existing amplitude scintillation prediction models often perform less satisfactorily when deployed outside the regions where they were formulated. This necessitates the need to evaluate the performance of scintillation models developed in one region using data data from other regions while documenting their relative errors. Due to its variation with elevation angle, frequency, other link parameters and meteorological factors, we employed three years (January 2016 to December 2018) of concurrently measured satellite radio beacons and tropospheric weather parameters to develop a location-based amplitude scintillation prediction model over the Earth-space path of Akure (7.17^oN, 5.18^oE), South-western Nigeria. The satellite beacon measurement used Tektronix Y400 NetTek Analyzer at 1 s integration time while meteorological parameters, namely; temperature, pressure and relative humidity were measured using Davis Vantage Vue weather station at 1 min integration time. Comparative study of the model’s performance with nine (9) existing scintillation prediction models indicates that the best and worst performing models, in terms of root mean square error (RMSE), are the Statistical Temperature and Refractivity (STN) and direct physical and statistical prediction (DPSP) models with values 11.48 and 51.03 respectively. Also, worst month analysis indicates that April, with respective enhancement and fade values of 0.88 and 0.90 dB for 0.01% exceedance, is the overall worst calendar month for amplitude scintillation.     

蜂窝状薄壁管拉伸时的内陷分析

蜂窝材料因具备广阔的结构设计空间,已广泛应用于航空航天等工业领域。探索蜂窝元胞的组织结构、组分性质与材料整体性能之间的关系有助于新型功能性材料的设计开发。本文从非平面Vertex模型的势能形式出发,结合蜂窝薄壁圆管拉伸时的轴对称特征,通过变分法得到了圆管拉伸时母线满足的控制方程,揭示了边界效应是蜂窝状薄壁圆管受拉时产生内陷的原因,并结合母线控制方程的若干特解,考察了非平面Vertex模型势能形式中包含的材料性质参数与圆管拉伸后内陷程度及范围的关系。研究还指出,蜂窝材料的材料性质包含在非平面Vertex模型的势能形式中,为蜂窝状材料整体性能的研究提供了一种新的思路。     

Solar sailing in realistic mode,based on a locally-optimal control laws

Ballistic design of solar sailing missions in the solar system is composed of defining the design parameters, the control programs, and the trajectories that provide performance goals of a flight. The use of a solar sail spacecraft imposes specific restrictions on mission parameters that include the degradation limit on the flight duration, the maximum temperature of solar sail's surface, the minimum distance from the Sun, the maximum angular velocity of the spacecraft's rotation and others.Many authors considered the impact of these restrictions on the design of the mission separately, but they used a sophisticated method of finding the exact optimal motion control or applied the most straightforward laws of motion control. This paper uses local-optimal control laws at the complete mathematical models of motion and functioning of solar sail spacecraft to describe a technique of designing interplanetary missions. The described method avoids the need to obtain an accurate optimal solution to the control problem and does not cause significant computational difficulties.     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

Initiation characteristics of oblique detonation waves from a finite wedge under argon dilution

In this paper, the flow field characteristics of Oblique Detonation Waves (ODWs) induced by a finite wedge under argon dilution are studied by solving the Euler equations with a detailed chemical model of hydrogen and air. First, the effects of the expansion waves, argon concentration, geometric parameters, and Mach number on the ODW are discussed. The results show that the changes of these parameters may make the oblique detonation not be initiated. Then, the ODW initiation criterion of the finite wedge is summarized, as the characteristic length of the induction zone L_C and the characteristic length of the oblique wedge L_W meet the condition L_C/L_W < 1, the initiation of the ODW occurs; otherwise, it does not occur. What’s more, the Constant Volume Combustion (CVC) theory is applied to study the characteristic length of induction zone. It is found that CVC theory is more suitable for the “smooth transition” type of ODW flow field, the theoretical and numerical characteristic length in induction regions are in good agreement. This work is of great significance for the design of oblique detonation engines for hypersonic vehicles.